Thermodynamic Properties of the S=l/2 Heisenberg Chain with Staggered 

Dzyaloshinsky-Moriya Interaction 
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Thermodynamic properties of the S=l/2 Heisenberg chain in transverse staggered magnetic field 
Hs and uniform magnetic field H x perpendicular to the staggered field is studied by the finite- 
temperature density-matrix renormalization-group method. The uniform and staggered magnetiza- 
tion and specific heat are calculated from zero temperature to high temperatures up to T/J — 4 
under various strength of magnetic fields from H% / 'J,H X / J — to 2.4. The specific heat and 
magnetization of the effective Hamiltonian of the Yb4As3 are also presented, and field induced gap 
formation and diverging magnetic susceptibility at low temperature are shown. 



I. INTRODUCTION 

The one-dimensional spin systems exhibit many in- 
teresting physical phenomena. Although these systems 
have only spin degrees of freedom, their ground state 
and thermodynamic behavior have diversity with vari- 
ous types of magnetic interactions. Recently, microscopic 
analysis JD on Yb4As3, a typical system showing field in- 
duced gap opening at low temperature, confirmed that 
its effective Hamiltonian is a S=l/2 Heisenberg model 
with Dzyaloshinsky-Moriya (DM) interaction. This mi- 
croscopic analysis supports previous scenario for the field 
induced gap opening due to the DM interaction [§|||Q, 
and stimulates more precise and detailed comparison be- 
tween experiment and theory. 

The Yb4As3 is a 4/ electron system, which undergoes 
charge ordering at 290K M . Below this temperature one 
of the four Yb ions becomes trivalent and forms one- 
dimensional (ID) chain along [111] direction [gQ. The 
Yb 3+ ion has one hole in the 4/ closed shell. The J = 7/2 
ground multiplet splits into four doublets by the crys- 
talline field effect. Thus the low temperature dynam- 
ics is described by an effective S=l/2 spin chain. The 
neutron scattering experiments on Yb4As3 actually con- 
firmed that the excitation spectrum is well described by 
the ID S=l/2 isotropic Heisenberg model |||j. 

Unusual features which can not be understood from 
the ID isotropic Heisenberg model are observed under 
applied magnetic fields. One is the upturn of the mag- 
netic susceptibility at low temperature, and the another 
is the formation of the excitation gap |l(j,[ll| . 

Recently, it was shown that these unusual features are 
caused by the staggered DM interaction, which generates 
effective staggered field under applied magnetic fields || . 
Since the staggered spin correlation length £(g=7r) of the 
Heisenberg model diverges at T — 0, drastic change oc- 
curs in the ground state. By using a mapping on to 
the sine-Gordon model, it has been shown that the ex- 
citation gap is formed by any finite staggered field pl[| . 
The analytic expressions of the gap, magnetic suscepti- 
bility, and specific heat of Yb4As3 are calculated within 
the sine-Gordon theory, and they are consistent with the 
experiments at low temperatures 0J . 

However the sine- Gordon model describes only the low 
energy physics of the effective one-dimensional S=l/2 
model. In the present paper, we directly calculate ther- 
modynamic quantities of the effective Hamiltonian of 
Yb4As3 over a wide range of temperatures and mag- 



II. MODEL AND METHOD 

The microscopic study on the Yb4As3 shows that 
its low temperature effective Hamiltonian is the S=l/2 
anisotropic Heisenberg model with the DM interaction: 



H = JJ2 {S?Si+i + cos2e(SfS? +1 + SfSf +1 )} 

i 

+ Jsm2eJ2(-^y(^x^+i) z +9±HY / Sf. (1) 

i i 

Here, H is the uniform external magnetic field perpen- 
dicular to the ID spin-chain. The DM interaction is elim- 
inated by rotating the spins in the staggered way in the 
x-y plane; 



Sf = 

Sf = -{-I) 1 sin^S 

QZ QZ 

Oa — Oa • 



cos e sf + {-iy sme sf, 



cose sj 



(2) 



After this transformation, the Hamiltonian is mapped on 
to 



H = J2_^ s» • Si 



■flx^E^ + flx^E^ 1 )^ 



(3) 



where H x — HcosO and H v s = Hsin9. This is equiv- 
alent to the isotropic Heisenberg model when external 
magnetic field is absent. Under finite external magnetic 
field H, both the effective uniform field H x and trans- 
verse staggered field H]j are applied to the Heisenberg 
spin-chain. 

In order 

to calculate thermodynamic quantities, we employ the 
finite-temperature density-matrix renormalization-group 
(finite-T DMRG) method g|@^|. Since the finite-T 
DMRG method iteratively expands the quantum trans- 
fer matrix in /3-direction restricting the number of basis 
states, the thermodynamic quantities are obtained for 
desired temperature. Thermodynamic quantities are ob- 
tained from the maximum eigenvalue and its eigenvector 
of the quantum transfer matrix, and the extrapolation on 
the system size is not needed [jfj). The numerical error 
is estimated from the eigenvalues of the density matrix 
which are truncated off, and thus the accuracy is sys- 
tematically improved by increasing the number of basis 



with the periodic boundary conditions. The thermody- 
namic limit is taken by using finite size scaling. 
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FIG. 1. Excitation gap of the ID Heisenberg model in the 
staggered field H^ . The uniform field H x is set to be zero. 
J = l.gj. = l. 
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FIG. 2. (a) Transverse staggered magnetization, (b) Uni- 
form magnetization in the transverse staggered field HH and 
the uniform field H x at T = 0. g± = 1. 
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FIG. 3. Transverse staggered magnetization at T 

i = l. 
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III. EXCITATION GAP AND MAGNETIZATION 

In Fig. 1 we show the excitation gap A as a function 
of staggered field H]j . The gap follows power law in the 
weak staggered field, and is well fitted by the exponent 
0.62 between E% ~ 0.002J and 2J. The power law de- 
pendence is also obtained in the field theory, which pre- 
dicts A = 1.78(.ffJ') 2/3 (-ln#]') 1/6 including logarithmic 
correction ||||] . This result is plotted in the figure by the 
dashed line. We find good agreement below H v a ~ 0.1J. 
In the region of intermediate strength of the staggered 
field Hj> ~ J, the gap slightly deviate from the pre- 
dicted power law, but still is quite close to the value 
A = 1.9(Hy)°- 62 . For strong field H% > 2 J, the spins are 
almost completely aligned to the staggered field and the 
gap is given by A = H^ + J as shown in the dotted line. 

The magnetization induced by the magnetic field is 
shown in Fig. 2. The transverse staggered magnetization 
M]{ at T — is shown in Fig. 2 (a), and the uniform 
magnetization M x is shown in Fig. 2 (b). For small stag- 
gered fields, the staggered magnetization monotonically 
increases with increasing uniform field until H x ~ 2 J. 
This result shows that the uniform field enhances the 
transverse staggered susceptibility. At H x ~ 2J, how- 
ever, the staggered magnetization sharply decreases. The 
decrease in MJ is related to the sharp increase in uniform 
magnetization M x as shown in Fig. 2 (b). The substan- 
tial change in the ground state at H x = 2 J is due to the 
singular behavior at the saturation of M x in the limit 
of weak Mf. For large H]j, both the staggered magne- 
tization and the uniform magnetization show monotonic 
change and no singular behavior is observed at H y ~ 2 J. 
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FIG. 4. Temperature dependence of transverse staggered 
magnetization at H]j = 0.01J. gi_ = 1. The excitation gap A 
is 0.11J at H x = 0. 
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FIG. 5. Specific heat in the staggered field H%. H x — 0. 
g± — 1. The excitation gap A is shown on the horizontal axis. 

In Fig. 3 we show the staggered field dependence of the 
staggered magnetization at fixed uniform field. We plot 

dl ^ M h o gH y 

dlogH! [ ° gUs 

to eliminate logarithmic dependence on the staggered 
field. In this figure we find linear dependence for weak 
staggered field. This result shows that M^ depends on 



H* as M» (x (Hv) a (logHV)P. The power law exponent 
a is determined by the slope in the figure, and a = 0.32 
is obtained at H x = 0. This is almost consistent with 
the result 1/3 obtained by the field theory @||. With 
increasing uniform field up to H x ~ 1.8 J, the power law 
exponent a decreases down to 0.24. The field theory also 
predicts the lowering of the power law exponent gj§] ■ 

The temperature dependence of the staggered magne- 
tization at Hy = 0.01J is shown in Fig. 4. At high 
temperature, M% is almost the same for H x < 1.6J. 
With decreasing temperature, MJ increases monotoni- 
cally and becomes almost constant below T/J ~ 0.2A. 
For H x = 2 J, the staggered magnetization is small com- 
pared with that for small uniform fields H x < 1.6J, but 
it continues to increase even at low temperatures below 
T/J ~ 0.01 due to the strongly reduced gap at H x = 2 J. 
For H x > 2 J, the spins are almost completely aligned to 
the uniform magnetic field and the temperature depen- 
dence is small. 
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FIG. 6. Specific heat in the transverse staggered field HJ* 
and the uniform field H x . g± = 1. 



IV. SPECIFIC HEAT 

Fig. 5 shows the specific heat at several strength of 
the staggered field H% under zero uniform field H x — 0. 
Since the excitation gap exists for finite H]j, deviation 
occurs from the T-linear dependence at T ~ A. The 
specific heat has a broad peak below T ~ A, and then 
exponentially decreases down to zero. As shown In Fig. 5 
(b), which show the specific heat divided by temperature, 
the sharp decrease occurs around the temperature T ~ 
0.2A with a broad peak at higher temperature T ~ 0.4A. 
The exponential dependence due to the gap is observed 
only below T - 0.2A. 

The effects of uniform field H x under a finite staggered 
field H]j = 0.0 U is shown in Fig. 6. The low temper- 
ature behavior below T ~ 0.2A is almost the same for 
H x < 1.6J. Since the low temperature behavior is char- 
acterized by the gap, this result shows almost constant 
gap for H x < 1.6 J. In contrast to the behavior below 



T ~ 0.2A, the broad peak at intermediate temperature 
strongly depend on the uniform field. The maximum of 
the specific heat around T/J ~ 0.4 at H x = becomes 
smaller and shifts to lower temperature with increasing 
H x . At H x = 1.6 the small dip appears at T/J ~ 0.2, 
and the main peak is divided into two peaks. With in- 
creasing H x , the peak at lower temperature shifts toward 
T = 0. The low temperature behavior is strongly mod- 
ified at H x — 2 J, and the peak at low temperature dis- 
appears for H x > 2J. With further increasing H x , the 
gap turns to increase with clear exponential temperature 
dependence at low temperature. 



V. APPLICATION TO YB4AS3 

The magnetization of Yb4As3 is calculated from the 
effective Hamiltonian Eq. (3) with rotating the spins 
back to the original coordinates. Thus the magnetiza- 
tion perpendicular to the ID chain of Yb4As3 is given 
by M = M x cos9 + Mfsin0. The calculated M are 
shown in Fig. 7 for various tan#, which is related to 
the ratio between the DM interaction (Dum) and J as 
sin 29 = Dom/J- The magnetization sharply increases 
near H = 0, and the exponent is 0.32 in the limit of 
H = 0. The magnetization curves gradually change its 
shape with decreasing tan#. By fitting the present re- 
sults with the experimentally obtained magnetization us- 
ing the experimentally estimated value of g± = 1.3, tan# 
is estimated to be 0.19 p7| . 

Temperature dependence of the magnetization divided 
by the magnetic field H is presented in Fig. 8. In this cal- 
culation we use the parameters tan 9 — 0.19, J = 2QK , 
and g± = 1.3 which reasonably reproduce susceptibility 
data M|. We also plot the susceptibility of ID Heisen- 
berg antiferromagnet (1D-HAF) for comparison. At high 
temperatures T > J the susceptibility is almost the same 
to that of 1D-HAF. With decreasing temperature, how- 
ever, the susceptibility of Yb4As3 largely increases espe- 
cially for small magnetic field due to the component of the 
transverse staggered susceptibility of the effective Heisen- 
berg model J4|. The susceptibility diverges at T = in 
the limit of H = 0. This is caused by the fact that 
the staggered correlation length of the Heisenberg model 
£(<7=7r) diverges at T = 0. 

The specific heat of Yb4As3 under the magnetic field 
perpendicular to the ID chain is presented in Fig. 9. The 
temperature dependence is similar to the results in Fig. 5, 
but in the present case, uniform field about four times 
larger than the effective staggered field induced by the 
DM interaction is present. Similarly to the results in 
Fig. 5 the specific heat C/T has a maximum around T = 
0.4A ~ 0.6A and exponentially decreases down to zero. 

It is highly desired that the theoretical results pre- 
sented in this paper will be compared with experimental 
results for a monodomain sample under magnetic fields 
perpendicular to the chain direction. 




H/J 
FIG. 7. Magnetization for various tan# at T — 0. g± — 1. 
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FIG. 8. Temperature dependence of the susceptibility of 
the effective Hamiltonian of YD4AS3. The magnetic field is ap- 
plied perpendicular to the ID chain. tan# = 0.19, J = 26i(', 
g±_ — 1.3. 1D-HAF is the susceptibility of ID Heisenberg 
antiferromagnetic chain in the limit of H — > 0. 




FIG. 9. Specific heat of the effective Hamiltonian of 
Yb4As3 in magnetic field perpendicular to the ID chain, 
tanfl = 0.19, J = 2QK, g± = 1.3. 
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